Biochimica et Biophysica Acta (BBA) - Biomembranes

1A-adrenoceptor (1A-AR) is a class A G-protein coupled receptor (GPCR) that stimulates smooth muscle contraction in response to adrenaline and noradrenaline. GPCRs exist in a dynamic equilibrium between multiple conformational states. Ligand binding induces structural rearrangements via conserved microswitches, which are thought to shift the equilibrium and trigger signalling. For structural and biochemical studies, GPCRs must be solubilised from the membrane, typically using detergent micelles. However, detergents can disrupt native dynamics of membrane proteins, potentially confounding experimental results. To address this, phospholipid bilayer mimetics such as nanodiscs and saposin nanoparticles (SNPs) have been developed to provide a more native-like environment. Thermostabilised 1A-AR serves as a GPCR prototype and can be expressed and isotopically labelled for NMR purposes. To investigate how membrane mimetics influence the conformational diversity of 1A-AR, we compared 1H 13C3-HMQC NMR experiments of 13CH3-Met labelled 1A-AR incorporated into either DDM, LMNG, or SNPs, in presence of ligands with varying efficacies. Several methionine residues are positioned near key microswitches, including M2035.57, located closed to the G protein binding site. Its resonance has been proposed as a readout of receptor conformational state, shifting with ligand efficacy. Spectra of 13CH3-Met labelled 1A-AR in LMNG closely resembled those in DDM with some temperature-dependent dynamic variation. In contrast, incorporation into SNPs led to a complete loss of M2035.57 signal, consistent with an intermediate exchange rate. These findings demonstrate that the membrane environment can profoundly influence conformational dynamics in GPCR NMR studies. Our results highlight the need to carefully consider membrane environment when interpreting NMR data and underscore the value of benchmarking against biologically relevant controls.

Liposomes are widely used as model membranes and nanoscale drug delivery systems, where cholesterol plays a key role in regulating bilayer structure and dynamics. However, how cholesterol concentration influences the structure and dynamics of liposome and how this influence is dependent on membrane curvature are not fully understood at the molecular level. In this work, coarse-grained molecular dynamics simulations using the MARTINI force field were employed to examine the concentration-dependent behavior of cholesterol in planar and curved membranes composed of cholesterol and unsaturated phospholipids, namely DOPC. More specifically, a planar lipid bilayer and an approximately 50-nm liposome were simulated to represent two extreme limits of small and large curvature, respectively. Increasing cholesterol concentration led to thicker membranes and reduced solvent exposure, consistent with cholesterols condensing effect. Membrane curvature enhanced interleaflet coupling and increased tail interdigitation relative to planar systems. Notably, DOPC flip-flop rate in spherical bilayers exhibited a non-monotonic dependence on cholesterol content, reflecting a balance between curvature-induced packing stress and cholesterol-driven ordering. These findings provide molecular-level insight into how cholesterol and curvature together shape the structure and dynamics of unsaturated lipid bilayers.

Nanodiscs and bicelles are widely used as membrane mimetics for structural studies of membrane-associated systems. Studies have reported that their magnetic alignment behavior and phase stability are highly sensitive to composition and temperature. In this study, we systematically investigate the effects of cholesterol on bicelle formation and magnetic alignment in DMPC + 0.2 glycyrrhizic acid (GA) systems using a combined 31P and 14N solid-state NMR experimental and simulation-based approach. Temperature dependent 31P NMR spectra reveal a clear transition from vesicle dominant to aligned bicelles/nanodsics phase, while 1N quadrupolar splitting and lineshape analysis provides quantitative insights into heterogeneous lipid bilayer populations, distinguishing large aligned nanodiscs (B(L)), small nanodiscs (B(S)), and isotropic/random components (B(R)). A strong correlation is observed between the 31P derived bicelle fraction and the 14N B(L) population, confirming that macroscopic alignment in the presence of an external magnetic field directly reflects the growth of large, well-ordered nanodiscs. Cholesterol is found to play a critical dual role by modulating membrane order and curvature. At low cholesterol concentration (0 to 5 mole percent), nanodiscs alignment occurs gradually with increasing temperature, while at higher cholesterol concentration (15 to 25 mole percent), the alignment is delayed and accompanied by broader spectral features, indicating structural heterogeneity. Notably, 10 mole percent cholesterol consistently provides the optimal balance, enabling efficient temperature dependent conversion to aligned bicelles while maintaining high B(L) populations (about 70-80 percentage) and minimal isotropic fractions. In contrast, higher cholesterol maintains significant B(S) and B(R) populations, even at elevated temperature. The 14N quadrupolar coupling (Cq is approximately 8.5 to 9.2 kHz for aligned nanodiscs) remains nearly invariant across compositions, showing that cholesterol does not change local headgroup dynamics but instead redistributes lipid populations. These findings establish a combined 31P and 14N solid -state NMR approach provides a valuable platform for quantitatively correlating membrane structure, dynamics, and alignment, offering practical guidelines for optimizing bicelle systems for high resolution solid-state NMR studies of membrane associated biomolecules.

The rising incidence of invasive fungal diseases emphasizes the need for novel therapeutic strategies, including membrane-targeting antifungal agents, which require representative lipid models for detailed molecular-level studies. In this work, we propose a consensus quinary fungal plasma membrane model based on lipidomic literature data, specifically PC:PE:PI:PA:PS phospholipid model with ratio of 44:29:13:8:6. Using a bottom-up approach, we characterized the biophysical properties of this system - with particular emphasis on mechanical parameters such as bending rigidity and area compressibility - by combining molecular dynamics simulations with experimental flicker-noise and ATR-FTIR spectroscopies. Furthermore, we investigated the effect of two key non-phospholipid components: ergosterol and triacylglycerols. Biophysical analysis revealed that DPPI and its specific interactions with DSPS induced the most substantial deviations in baseline membrane parameters, particularly area per lipid, membrane thickness, and area compressibility, while DSPS influenced bending rigidity change and DLiPA primarily affected lipid packing defects. In addition, ergosterol and TGs were found to influence all of the investigated parameters to different degree. Notably, the overall biophysical profile of the proposed FPMM closely mimicked that of natural vesicles derived from yeast lipid extracts, establishing this model may provide a reliable platform for studying fungal membrane biophysics and lipid-targeting interactions.

Magnetic treatment in tomato seeds (Solanum lycopersicum L.) has been studied as a biotechnological technique to induce a reduction in germination times and enhance plant development. However, the modified cellular mechanisms involved in the reduction of germination times or the improvement of development parameters are not yet clearly established. To explore a possible altered cellular mechanism, the effect of homogeneous static magnetic fields on the structure of the cyclic nucleotide-gated channel 6 (CNGC6), the modification in the organization of POPC lipids in the plasma membrane, and changes in calcium ion mobility were evaluated. For this purpose, coarse-grained molecular dynamics simulations were performed using the Martini 3 model in GROMACS, applying five different magnetic flux densities (0.000, 0.001, 0.010, 0.100, 1.000, and 10.000) T over 1 000 ns. The results showed an anisotropic effect in the longitudinal direction of the protein, which generated heterogeneous behavior among the chains of the homotetramer; this altered the conformation of the CNGC6 channel and modified the pore bottleneck. In contrast, no significant changes were observed in the conformational order of the POPC phospholipid chains. As a preliminary, single-replicate exploratory study, these results suggest that homogeneous static magnetic fields may induce specific structural modifications in the CNGC6 ion channel of Solanum lycopersicum L. without compromising the integrity of the lipid bilayer or the dynamics of ion transport within the analyzed timescale; these preliminary findings provide a molecular-level structural basis for future experimental and computational investigations of magnetic field effects on plant cyclic nucleotide-gated channels.

Neutron scattering is a powerful technique to investigate atomic structures and molecular dynamics of proteins at the nano-scale. When it comes to dynamics, incoherent and coherent scattering respectively provide information on the single and collective dynamics of nuclei. In proteins, hydrogen has the highest incoherent cross-section, and it is common practice to overlook the contribution of coherent terms stemming from all nuclei. However, the fast collective dynamics of heavier nuclei could also be studied if coherent scattering and incoherent scattering were experimentally separated. The recent advent of polarized neutron spectroscopy with sufficient flux and energy resolution has made it possible, and opens new perspectives to investigate the relative importance of coherent scattering and the information it provides on biological samples. The present study reports on the use of polarized quasi-elastic neutron scattering (QENS) and the application of a minimalistic model adapted to both individual and collective dynamics. Using a perdeuterated green fluorescent protein as a model globular protein, the study provides an interpretation of the dynamical parameters obtained with QENS, and a comparative study of the Elastic Coherent and Incoherent Scattering Factor. Based on both experiments and calculations, we discuss the relative importance of distinct and self components of coherent scattering, which is often wrongly assumed to be representative of collective dynamics only. The results highlight the current impediments rendering complicated a straightforward analysis of fast collective dynamics in hydrated protein samples.

Neuroendocrine cells communicate with other cells by releasing neurotransmitters or hormones by exocytosis, which involves SNARE-mediated fusion between secretory vesicles and the plasma membranes of the secreting cells. In neurons two plasma membrane SNARE proteins, Syntaxin-1a and SNAP25, join with the vesicle membrane SNARE protein Synaptobrevin-2 to form a four-helix bundle, which drives membrane fusion. The assembly of these SNAREs, which is highly orchestrated in cells, has been intensely studied in solution using fragments of the SNARE proteins without their transmembrane domains or lipid anchors. However, in cell and model membranes, Syntaxin and SNAP25 are known to oligomerize and cluster, and little is known about how clustering affects their incorporation into SNARE complexes. In cells, the SM protein Munc18 has been implicated in aiding secretory vesicle docking and facilitating SNARE complex assembly through its interactions with Syntaxin. To understand how Munc18 orchestrates SNARE complex assembly on membranes, we employed protein reconstitution in model membranes as well as biochemical and biophysical assays to show that lipid-dependent oligomerization of Syntaxin affects Munc18-Syntaxin binding and SNAP25 insertion into the plasma membrane acceptor SNARE complex. We showcase the consequences of the different modes of Munc18-Syntaxin and SNAP25 interaction on Syntaxins oligomerization and orientation relative to the membrane surface, as well as on docking and fusion of purified insulin granules. We also determined low-resolution structures by cryoEM in nanodiscs and on the surface of proteoliposomes of membrane-bound assembly states of Munc18/Syntaxin and Munc18/Syntaxin/SNAP25 complexes.

The Epstein-Barr virus (EBV) is a highly prevalent virus worldwide that is associated with several lymphoid and epithelial malignancies. However, extensive research on EBV integral membrane proteins BILF1, LMP1 and LMP2, has been scarce due to their hydrophobic transmembrane domains. Our study applies the QTY code (glutamine, threonine, tyrosine) to design water-soluble analogs of BILF1, LMP1 and LMP2 with reduced hydrophobicity, where we systematically replaced hydrophobic amino acid residues leucine (L), isoleucine (I), valine (V), and phenylalanine (F) with structurally similar polar residues glutamine (Q), threonine (T), and tyrosine (Y). We retrieved their native sequences from UniProt, identified transmembrane domains using Protter, then performed QTY design through the Protein Solubilizing Server (PSS). We then predicted native and QTY structures using in silico prediction tools AlphaFold3, ColabFold, and Boltz-2. Our analyses demonstrate that despite significant protein sequence replacements in their transmembrane domains (54.15%-61.59%) and increased intrinsic solubility, the QTY analogs exhibited minimal changes in isoelectric point (0.00-0.15 decrease) and molecular weight (0.7-1.2 kDa increase). Additionally, structural superpositions between QTY analogs and native structures using PyMOL yield low RMSD values (0.217[A] -1.202[A]). Our results demonstrate the QTY codes ability to design detergent-free analogs of BILF1, LMP1 and LMP2 with substantially reduced hydrophobicity and aggregation propensity whilst preserving native-like structures. Our results may facilitate protein characterization studies, therapeutic research on EBV, and other protocols that typically require protein solubilization.

Cowpox virus CPXV012 inhibits MHC class I antigen presentation by interfering with TAP-dependent peptide transport, but its membrane-dependent structural organization and dynamic behavior remain incompletely defined. Here, we investigated the conformational properties of CPXV012 in membrane-mimicking environments and in a model of the CPXV012-TAP complex. CPXV012 was divided into three peptide constructs corresponding to the N-terminal cytosolic region, transmembrane segment, and C-terminal ER-luminal domain. The peptides were analyzed by circular dichroism spectroscopy, multidimensional NMR spectroscopy, and molecular dynamics simulations, and the resulting structural information was integrated into a full-length CPXV012 model. CD spectra showed that CPX-E1 and CPX-C2 are predominantly disordered in aqueous solution but acquire ordered, mainly -helical features in DPC micelles. NMR analysis in DPC-d38 micelles provided residue-level assignments and structural restraints supporting restrained structure calculations for both peptides. In three independent 1 {micro}s molecular dynamics simulations of the CPXV012-TAP complex, CPXV012 preserved a reproducible two-helical organization. The N-terminal/transmembrane region behaved as a relatively stable structural element, whereas the ER-luminal segment showed greater local flexibility. Interface analysis indicated that CPXV012 contacts both TAP1 and TAP2, with recurrent interactions concentrated in the luminal Y47-I69 region and involving polar and charge-complementary contacts. These results support a model in which membrane-associated structuring positions CPXV012 for TAP recognition, while the flexible ER-luminal region forms the main TAP-interacting surface. This structural framework complements existing functional models of CPXV012-mediated TAP inhibition.

The BCL-2 family proteins are key regulators of apoptosis, functionally divided in pro- and anti-apoptotic proteins, with a third group acting as regulators. Their ability to partition between the cytosol and intra-cellular membranes (essentially the mitochondrial outer membrane) is a primary regulator of their functions. A second contributor is their ability to form homotypic complexes (pro-pro or anti-anti) or heterotypic complexes (pro-anti). If the structures of monomeric cytosolic members have largely been characterized, the functional and structural study of membrane-embedded proteins remains incomplete. Unlocking this knowledge is expected to enable evaluating new therapeutic strategies to either activate pro-apoptotic members, or inactivate anti-apoptotic ones. Lipid bilayer nanodiscs and improved cell-free protein synthesis have provided the technical breakthrough to achieve the description at the atomic level of conformations and higher order assemblies of these proteins in their membrane-associated states. Here we describe detailed and straightforward protocols for generating nanodisc-inserted members of the Bcl-2 family, through the example of anti-apoptotic Bcl-xL, and pro-apoptotic Bax and Bak. Full-length, untagged proteins are expressed from bacterial extracts in the presence of pre-assembled nanodiscs to allow co/post-translational insertion in lipid bilayer, followed by affinity chromatography purification. A more detailed characterization is presented for Bak, to exemplify structural and mechanistic studies enabled by these methods. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=139 SRC="FIGDIR/small/745005v1_ufig1.gif" ALT="Figure 1"> View larger version (35K): org.highwire.dtl.DTLVardef@5da1d1org.highwire.dtl.DTLVardef@12aca96org.highwire.dtl.DTLVardef@5a3e73org.highwire.dtl.DTLVardef@ba009d_HPS_FORMAT_FIGEXP M_FIG C_FIG

We present a comparative study of the inner membrane of three Gram-positive bacterial strains, namely S. aureus, S. epidermidis and N. lacusekhoensis. A lipidomics study is used to obtain the lipid architecture and composition for S. epidermidis found in the skin microbiome and N. lacusekhoensis, an extremophile present in halophilic and alkophilic environments. Differences between the strains arise from both the lipid architecture and the cardiolipin content varying from 5% in S. aureus to 85% in N. lacusekhoensis. We develop coarse grained (CG) Martini-3 membrane models which reproduce structural properties such as membrane area, thickness, density distributions as well as ion-correlations with all-atom models. Inter-lipid correlations reveal a homogeneous distribution of lipids in the membranes despite the wide variation in lipid types and composition. Mechanical properties such as the area stretch modulus increased with cardiolipin content, however the bending modulus has a more complex dependence on membrane charge and lipid type. Using the CG models we evaluate the insertion free energies for four widely used antimicrobial molecules. Entry barriers for thymol and methylparaben arise from the charge density modulation at the membrane headgroups due to counterion condensation. The entry mechanisms of the antimicrobial peptide cecropin-melittin-15 (CM15) and the preservative molecule ethyl-lauroyl-arginate (ELAR) are found to be similar across all three strains. We also illustrate the manner in which the extremophilic strain, N. lacusekhoensis with its high cardiolipin content, modulates the partitioning kinetics of the antimicrobial molecule thymol with pH and salt. Our study reveals that membrane properties are largely conserved across the three model membranes. The molecular models and insights emerging from the present work should aid in the development of novel antimicrobials against Gram-positive strains.

The perilipins generally represent the most abundant lipid droplet (LD) surface proteins in mammalian cells and can target LD subpopulations within the same cell. They are characterized by a conserved organization of disordered and folded regions, as well as a number of divergent features, which contribute to differences in perilipin function and LD targeting. Here, we focus on the C-terminal 4-helix bundle (4HB) domain that is present in all perilipins except for PLIN1. Using biochemical and in silico approaches, we show that the 4HB of PLIN3 is a stably folded domain and interacts with lipid surfaces in vitro and with LDs in model cells. The {beta}-subdomain at the bottom of the helical bundle is required for the binding to LDs, but not for the 4HB stability, suggesting that this region may promote direct interaction with the LD surface. In agreement, the 4HB of PLIN4, which does not contain an {beta}- subdomain, does not bind to LDs. Overall, our work shows that small differences in perilipin structural features impact their differential targeting to LDs.

Abstract textMembrane protein function depends critically on the surrounding lipid bilayer, which supports specific lipid-protein interactions and imposes constraints through its ensemble properties. Ion-transporting ATPases undergo large conformational changes that perturb these interactions, and the high chemical potential of water can further drive conformational transitions through changes in internal hydration. Here, we investigated the relationship between membrane environment, internal hydration, and catalytic activity in the P1B-type copper ATPase LpCopA from Legionella pneumophila. BADAN-labeled LpCopA mutants revealed distinct basal hydration states near the conserved canonical binding site (CBS) cysteines C382 and C384, with C382 residing in a more hydrated environment than C384. Using the osmolyte PEG-1500, the spectral response of BADAN indicated an average internal hydration volume of about 800 [A]3 for LpCopA in E. coli lipid-doped mixed micelles (MMs), corresponding to approximately 25 water molecules. In nanodiscs (NDs), less than half of this volume responded to osmotic pressure, consistent with compaction of the transmembrane helical bundle under membrane lateral pressure. Remarkably, ATPase activity measured in different lipid reconstitution systems scaled with the extent of internal hydration, with diisobutylene/maleic acid lipid particles (DIBMALPs) imposing the tightest transmembrane compaction and lowest basal activity. These data indicate that the packing density and internal hydration of the transmembrane domain of LpCopA are strongly modulated by the membrane environment. The energetic estimates further support a model in which hydration-dependent expansion of the transmembrane domain against membrane lateral pressure contributes to the free-energy barrier of ATP hydrolysis.

Reconstructing time-resolved inter-residue distance distributions during protein functional dynamics in the solution state is known to be a difficult and important problem. This article presents a technique for extracting spin-spin (as a proxy for residue-residue) distance distributions on doubly-spin-labeled proteins from rapid-scan time-resolved Gd-Gd electron paramagnetic resonance (rs-TiGGER) spectra recorded near room temperature in solution at 240 GHz. We use a best-fit technique that convolves a dipolar kernel matrix with an intrinsic, non-dipolar-broadened (single-labeled) spectrum. The kernel incorporates the effect of solution-state tumbling on the dipolar broadening using a correlation function that bridges the static and rapidly tumbling regimes. We apply the technique to AsLOV2, a protein domain with a dark-state crystal structure that is well-known from X-ray crystallography, but a less well-characterized and disordered tertiary structure that manifests after photoactivation at 450 nm. Informed by principal component analysis, we assume that the underlying distance distribution may be approximated by a sum of two Gaussian distributions. The fits returned time-resolved, light-activated populations with mean distances of [Formula] (dark) and [Formula] (lit) in the wild type, and [Formula] (dark) and [Formula] (lit) in an N414Q mutant, with nearly complete unfolding (within fit uncertainty) of the active, light-sensitive fraction. The extracted distance distributions and their accompanying uncertainties are consistent within uncertainty with molecular dynamics simulations of the equilibrated protein structure.

Production of recombinant proteins frequently yields inclusion bodies that must undergo refolding to yield active protein. Here, we optimized the refolding conditions for the recombinant leucyl aminopeptidase (rPepL) from Lactocaseibacillus casei expressed in inclusion bodies from E. coli. Several chemical additives were assessed for how well they facilitated an increase in refolding efficiency. The best, 0.5 M L-arginine, yielded 50.8% refolding. The addition of stabilizers, such as sucrose and glycerol, with L-arginine further increased yields to 85%. Urea at lower concentrations (0.25-0.5 M) also facilitated an increase in the refolding yield when co-added with L-arginine, whereas guanidinium chloride inhibited it. Sugars and polyols exhibited dose-dependent effects, with ranges for optima also defined. Fluorescence spectroscopy verified enhancements in the refolding under the optimized conditions. Molecular dynamics simulation under mixed solvent conditions provided atomic insights about stabilizing interactions that are likely to facilitate increased refolding. The results show that a series of aggregation suppressors and protein stabilizers can, in a collaborative way, increase the refolding efficiency for the recombinant proteins from the inclusion bodies. The protocol with the optimization using the additives L-arginine, sucrose, and glycerol is an efficient method for the production of active rPepL. This article outlines the best refolding method to recover recombinant leucyl aminopeptidase from inclusion bodies of E. coli using L-arginine combined with sucrose and glycerol. The combined experimental observations and computational simulations elucidate the molecular process of additive-induced stabilization, which elucidates how aggregation inhibition and hydrogen-bonded stabilization act synergistically. The results presented herein answer both mechanistic understanding and experimental guidance for improving protein refolding.

Small-angle X-ray or neutron scattering (SAXS/SANS) analysis of large unilamellar vesicles (LUVs) is often limited by high-dimensional bilayer models and the lack of dedicated, statistically rigorous workflows. Here, we introduce SAS_MoCa, an open-source Python package that integrates a compositional scattering density profile (SDP) description of lipid bilayers with a separated form factor (SFF) treatment of vesicle size and polydispersity, and couples these highly parameterized models to an adaptive thermodynamic simulated annealing algorithm formulated within a constrained Bayesian framework. SAS_MoCa enables users to incorporate quantitative prior information from, e.g., previous SAXS/SANS studies, dynamic light scattering, NMR, or molecular simulations, and returns full posterior parameter distributions, uncertainties (reported as medians and median absolute deviations) and correlations even from single SAXS curves. Validation on POPC, POPE and DMPC SAXS-only data demonstrates that the method yields reproducible structural parameters with uncertainties comparable to joint SAXS/contrast-variation SANS analyses. The modular architecture of SAS_MoCa facilitates extension to additional lipid systems and future joint SAXS/SANS or SANS-only applications.

Cardiolipin (CL) is a phospholipid found in the inner mitochondrial membrane (IMM) where it increases the efficiency of ATP regeneration. We have investigated the hypothesis that this increase may result in part from CL concentrating H+ at the IMM surface through electrostatic interactions as the CL polar head is a dianion at physiological pH. To this aim, we compared the concentrations and movements of H+ at the surface of giant planar phosphatidylcholine (PC) membranes and 20% CL enriched PC membranes by recording their surface pH with the membrane-grafted pH probe fluorescein DHPE. CL enrichment of the membranes increased their surface H+ activity by a ~4 factor. Moreover, we observed non-gaussian spatial H+ concentration profiles with distance from a point H+ source with both PC and CL membranes suggesting that both lipids also induce interactions between probe molecules. A whole bath pH variation revealed that these interactions allow the traveling of reversible acidification fronts with constant speed over the membrane between high and low pH states. A reaction-diffusion model of these observations suggests that membranes support these fronts through a mechanism of autocatalytic (de)protonation of the membrane surface. In mitochondria, these fronts would result in transitions between high and low pH states, the low one having a larger H+ concentration in CL-enriched regions of the IMM. Such an increase at the inner leaflet of the IMM may increase efficiency of the respiratory chain whereas the increase at the outer leaflet may boost the ATP synthase rate.

Lipid droplets (LDs) are cellular organelles responsible for lipid storage and metabolism. The mechanism of biogenesis of LDs involves phase separation of neutral lipids from the surrounding phospholipids, which generates oil lenses embedded in lipid bilayers, also known as nascent LDs. As nascent LDs grow, at some point they bud out of the bilayer, forming nearly spherical droplets. Nascent LDs have different propensity to bud, and it has been proposed that their shape provides information on such propensity; however, LD shape is difficult to determine experimentally. Here we studied the shape of lipid droplets using MD simulations at the coarse-grained level, and compared it to the predictions by an established theory. Our general system setup features an oil lens embedded into a flat, periodic bilayer. We found that the shape of simulated nascent LDs resembles a spherical cap (i.e., it has constant curvature over most of the surface), in excellent agreement with the theory, already for very small droplet sizes. The aspect ratio (height/radius) of nascent LDs increases with increasing LD volume, increasing membrane softness, and increasing surface tension between oil and water, also in agreement with theoretical predictions; however, it remains lower than 1 (i.e., the ratio for a sphere) for LDs of up to 40 nm in diameter. Fitting the simulated LD shapes with a theoretical shape equation suggests that a non-zero surface tension is present in both the monolayer and in the bilayer region. The existence of a relatively high surface tension in the bilayer region is confirmed by local stress calculations, and indicates that the periodic system setup does not reproduce the properties of nascent LDs in the endoplasmic reticulum, where the bilayer tension is two orders of magnitude lower. However, the simulations provide a microscopic view into the properties of droplet embedded vesicles.

Live cell plasma membranes show spatially heterogeneous liquid-ordered (Lo)-like and liquid-disordered (Ld)-like regions similar to the co-existing Lo/Ld phases observed in lipid vesicles. The Lo-like regions are relatively less hydrated and less polar due to tight packing of the membrane components compared to the Ld-like regions. The steady-state fluorescence spectra of Di-4-ANEPPDHQ (Di-4), a widely used polarity-sensitive probe, is blue or red shifted when solvated in less polar (Ld- like) or more polar (Lo-like) regions respectively. However, quantification of Di-4 fluorescence in blue and red channels for the evaluation of membrane phase state suffers from the lack of specific wavelength choice for these two channels and Di-4s relatively higher concentration in Ld phase (red channel) due to its partitioning preference. To address these issues, we employed fluorescence lifetime of Di-4, a concentration independent photophysical parameter, to understand membrane biophysical properties. The fluorescence lifetime of Di-4 in lipid vesicles exhibits Arrhenius-like temperature dependence. Centred around this energetic feature of Di-4 photophysics, we developed a novel analytical module, namely excited state relaxation activation energy (ESRAct), that serves as an intrinsic descriptor of the membrane nano-environment sensed by this probe. We show that the ESRAact value scales with increasing disorder in nanoscale phase separation (i.e., ESRAct of pure Ld > mixed Ld/Lo > pure Lo phase). We then extended its applications to giant plasma membrane vesicles (GPMVs) isolated from MCF-7 cells and found that these vesicles exhibit nanoscale Lo/Ld co-existing phase within 16-37{degrees}C. We envisage wide applications of ESRAct to delineate plasma membrane phase behavior as well as general photophysical studies on other newly designed polarity-sensitive probes.

Interactions between proteins and glycosphingolipids (GSLs) regulate various cellular processes and altered GSL metabolism contributes to numerous diseases. The diverse glycan headgroups and ceramide backbones of GSLs shape membrane organisation, fluidity, curvature, and tension. As protein recognition frequently depends on both glycan specificity and the organisation of GSLs within the membrane, these interactions remain challenging to characterise in vitro. Here, we introduce FLiPA (Fluorescent Liposome Plate Assay), a versatile method that utilises fluorescent agarose-embedded giant liposomes for the quantitative analysis of protein-GSL interactions. By enabling systematic control of membrane and buffer composition, FLiPA provides an accessible and robust platform for dissecting the molecular determinants of protein-GSL interactions, including the roles of cholesterol, membrane order, protein oligomerisation and ionic strength.